Predicting rogue waves in random oceanic sea states 
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^ I Using the inverse spectral theory of the nonlinear Schrodinger (NLS) equation we 

correlate the development of rogue waves in oceanic sea states characterized by 
the JONSWAP spectrum with the proximity to homoclinic solutions of the NLS 
equation. We find in numerical simulations of the NLS equation that rogue waves 
develop for JONSWAP initial data that is "near" NLS homoclinic data, while 
CN ' rogue waves do not occur for JONSWAP data that is "far" from NLS homoclinic 

^ ■ data. We show the nonlinear spectral decomposition provides a simple criterium 

CN ' for predicting the occurrence and strength of rogue waves (PACS: 92.10.Hm, 

^ ■ 47.20.Ky, 47.35+i). 

' Introduction Rogue waves are rare, large amplitude waves whose heights exceed 2.2 times 



O 



the significant wave height of the background sea. One of the proposed mechanisms for the 
development of rogue waves in deep water is nonlinear focusing due to the Benjamin- Feir 
! (BF) instability.^'^ The BF instability is a modulational instability in which a uniform train 

^ I of surface gravity waves is unstable to a weak amplitude perturbation. The BF instability 
I is described approximately by the focusing nonlinear Schrodinger (NLS) equation ^ 

. ; iut + u^x + '2\u\\ = 0, (1) 

and in the simplest setting homoclinic orbits of the unstable Stokes solution of the NLS 
equation have been used for modeling rogue waves. ^'^ Homoclinic solutions of the NLS 
equation, obtained when two or more unstable modes are present, can be phase modulated 
to provide striking examples of wave amplification where the amplification is due to both 
the BF instability and the additional phase modulation.^ 

The NLS equation is the leading order equation in a hierarchy of envelope equations 
and is derived from the full water wave equations under the assumption of a narrow 0{e) 
banded spectrum. This bandwidth constraint limits the applicability of the NLS equation 
in 2D as it results in energy leakage to high wave number modes. ^ The broader bandwidth 
NLS (BBNLS) equation, obtained by assuming the bandwidth is 0{y/e) and by retaining 
higher order terms in the asymptotic expansion for the surface wave displacement, has been 
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successful in reducing the energy leakage.^ An alternate approach is to "enhance" the NLS 
equation with exact hnear dispersion, whereby the equation has improved bandwidth res- 
olution and stability properties/ All of these higher order equations, whether narrow or 
broader bandwidth or "enhanced", may be viewed as perturbations of the NLS equation. 
Homoclinic orbits of the Stokes wave have been shown to persist for the BBNLS equation. ^'^ 
This persistence result suggests homoclinic solutions of the NLS equation may be significant 
in modeling rogue waves for random oceanic states. 

Onorato et al.^ examined the generation of extreme waves for typical random oceanic 
sea states characterized by the Joint North Sea Wave Project (JONSWAP) power spectrum. 
In numerical simulations of the NLS equation it was found that rogue waves occur more 
often for large values of the Phillips parameter a and the enhancement coefficient 7 in the 
JONSWAP spectrum. Even so, they observed that large values of a and 7 do not guarantee 
the development of extreme waves. 

In this paper we clarify the dependence of rogue wave events on the phases in the "ran- 
dom phase" reconstruction of the surface elevation (see eqn. 0). We find that the phase 
information is as important as the amplitude and peakedness of the wave (governed by a 
and 7) when determining the occurrence of rogue waves. Random oceanic sea states char- 
acterized by JONSWAP data are not small perturbations of Stokes wave solutions. As a 
consequence, it is difficult to investigate the generation of rogue waves in more realistic sea 
states using a linear stability analysis (as in the Benjamin- Feir instability). Our approach 
is based on the NLS equation and its inverse spectral theory, used to examine a nonlinear 
mode decomposition of JONSWAP type initial data. Such analysis allows us to determine 
the nonlinear mode content of the data and the proximity (measured in terms of a parameter 
5) to instabilities and homoclinic solutions of the NLS equation. 

Our main results are: 1) JONSWAP data can be quite near data for homoclinic orbits of 
complicated A^-phase solutions. For fixed values of a and 7 in the JONSWAP spectrum, as 
the phases in the initial data are randomly varied, the proximity 5 to homoclinic data varies. 
2) In several hundred simulations of the NLS, where the parameters and the phases in the 
JONSWAP initial data are varied, we find that rogue waves develop for JONSWAP data 
that is "near" NLS homoclinic data, while rogue waves do not occur for JONSWAP data 
that is "far" from NLS homoclinic data. Consequently, we find that the nonlinear spectral 
decomposition provides a simple criterium, in terms of the proximity to homoclinic solutions, 
for predicting the occurrence and strength of rogue waves. This is the first time homoclinic 
solutions have been correlated with rogue waves for realistic oceanic conditions. 
Random oceanic sea states To examine the generation of rogue waves in a random sea 
state, we note that the surface elevation rj is related to m, the solution of the NLS equation. 



where C„ is the amplitude of the nth component with wave number kn = {n — l)k, k = 2tt/L, 
and random phase 0„, uniformly distributed on the interval (0, 27r). The spectral amplitudes. 
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Cn = —i\/2Sn/ L, are obtained from the JONSWAP spectrum:^ 

/ f - f „ \ ^ 

(3) 
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Here /o is the dominant frequency, determined by the wind speed at a specified height above 
the sea surface, ctq = 0.07 (0.9) for / < /o (/ > /o) and fn = n/L is the wave frequency. 
The parameter 7 is the peak-shape parameter; as 7 is increased, the spectrum becomes 
narrower about the dominant peak. For 7 > 1 the wave spectra continues to evolve through 
nonhnear wave-wave interactions even for very long times and distances. It is in this sense 
that JONSWAP spectra describe developing sea states rather than a fully developed sea. The 
scale parameter a is related to the amplitude and energy content of the wavefield. Based 
on an "Ursell number", the ratio of the nonlinear and dispersive terms of the NLS equation 

in dimensional form, the NLS equation is considered to be applicable for 2 < 7 < 8.^ 
Typical values of alpha are .008 < a < .02. 

We examine a nonlinear spectral decomposition of the JONSWAP initial data, which 
takes into account the phase information This decomposition is based upon the inverse 
scattering theory of the NLS equation, a procedure for solving the initial value problem 
analogous to Fourier methods for linear problems. We find that we are able to predict the 
occurrence of rogue waves in terms of the proximity 6 to distinguished points of the discrete 
spectrum. We briefly recall elements of the nonlinear spectral theory of the NLS equation. 
Floquet spectral theory The integrability of the NLS equation ([Q) is related to the fol- 
lowing pair of linear systems (the so-called Lax pair) 



n"". r =0, £''v=o, (4) 




where D denotes the derivative with respect to x, A is the spectral parameter and (j) is 
the eigenfunction.^ These systems have a common nontrivial solution (f){x,t]X), provided 
the potential u{x,t) satisfies the NLS equation. C^^^ is not specified explicitly as it is not 
implemented in our analysis. 

The first step in solving the NLS using the inverse scattering theory is to determine 
the spectrum a{u) = |A G C | ^C^^V = 0^ \4>\ bounded Vxj of the associated linear operator 
C^^\ which is analogous to calculating the Fourier coefficients in Fourier theory. For pe- 
riodic boundary conditions, u{x + L,t) = u{x,t), the spectrum of u is expressed in terms 
of the transfer matrix M{x + L; u, A) across a period, where M{x; u, A) is a fundamen- 
tal solution matrix of the Lax pair Q. Introducing the Floquet discriminant A(u, A) = 
Trace [M(x + L; u, A)], one obtains'^ 

a{u) = {A G C I A(m, A) G R, -2 < A(m, A) < 2} . (5) 

The distinguished points of the periodic/antiperiodic spectrum, where A{X,u) = ±2, are: 
(a) simple points {A^" | dA/d\ ^ 0} and (b) double points {Aj | dA/d\ = 0, d^A/dX'^ ^ 0}. 
The Floquet discriminant functional A{u, A) is invariant under the NLS flow and encodes 
the infinite family of constants of motion of the NLS (parametrized by the Xj). 
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The Floquet spectrum (jSJ of a generic NLS potential consists of the entire real axis 
plus additional curves (called bands) of continuous spectrum which terminate at the simple 
points A^. A^-phase solutions are those with a finite number of bands of continuous spectrum. 
Double points arise when two simple points have coalesced and their location is important. 

Using the direct spectral transform, any initial condition or solution of the NLS can 
be represented in terms of a set of nonlinear modes. The spatial structure and dynamical 
stability of these modes is determined by the order and location of the corresponding Xj as 
follows: (a) Simple points correspond to stable active degrees of freedom, (b) Double points 
label all additional potentially active degrees of freedom. Real double points correspond to 
stable inactive (zero amplitude) modes. Complex double points are associated with all the 
unstable active modes and label the corresponding homoclinic orbits. 

Figure 1 shows the spectrum of a typical unstable A^-phase solution. There are N bands 




Figure 1: Spectrum of an unstable A^-phase solution. 

of spectrum determined by the 2N simple points A^. The 2M complex double points A^ 
indicate that the solution is unstable and that there is a homoclinic orbit. The simple 
periodic eigenvalues are labeled by circles and the double points are labeled by crosses. An 
example of a spectrum for a nearby semi-stable A^-phase solution where the complex double 
point is split 0{e) is given in Fig. 2(a). 

Explicit formulas for the A^-phase solutions, Q{6i, ...,6n), are obtained in terms of the 
simple spectrum. The phases evolve according to 6j = kjx + Qjt + 6j, kj = 27mj/L, where 
Kj and Qj are determined by A^ (since the spectrum is invariant and Qj are constants). 
For a given A^-phase solution, the isospectral set (all NLS solutions with the same spectrum) 
comprises an A^-dimensional torus characterized by the phases 6j. If the spectrum con- 
tains complex double points, then the A^-phase solution may be unstable. The instabilities 
correspond to orbits homoclinic to the A^-phase torus. 

JONSWAP data and the proximity to homoclinic solutions of the NLS In the 

numerical simulations the NLS equation is integrated using a pseudo-spectral scheme with 
256 Fourier modes in space and a fourth order Runge-Kutta discretization in time {At = 
10^'^). The nonlinear mode content of the data is numerically computed using the direct 
spectral transform described above, i.e. the system of ODEs (j3)) is numerically solved to 
obtain the discriminant A. The zeros of A it 2 are then determined with a root solver based 
on MuUer's method. The spectrum is computed with an accuracy of O{10~^), whereas the 
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spectral quantities we are interested in range from (9(10^^) to (9(10^^). 

Complex double points typically split under perturbation into two simple points, \±, thus 
opening a gap in the band of spectrum (see Fig. 2(a)). We denote the distance between these 
two simple points by 5(A+, A_) = |A+ — A_| and refer to it as the splitting distance. We use 
6 to measure the proximity in the spectral plane to homoclinic data, i.e. to complex double 
points and their corresponding instabilities. Since the NLS spectrum is symmetric with 
respect to the real axis and real double points correspond to inactive modes, in subsequent 
plots only the spectrum in the upper half complex A-plane will be displayed. 

We begin by determining the spectrum of JONSWAP initial data given by Q for various 
combinations of a = .008, .012, .016, .02, and 7 = 1,2,4,6,8. For each such pair (7,0;), 
we performed fifty simulations, each with a different set of randomly generated phases. As 
expected, the basic spectral configuration and the number of excited modes depended on the 
energy and the enhancement coefficient a and 7. However, the extent of the dependence of 
the spectrum upon the phases in the initial data was surprising. 

As a typical example of the results. Figs. 2(a) and 3(a) show the numerically computed 
nonlinear spectrum of JONSWAP initial data when 7 = 4 and a = .016 for two different 
realizations of the random phases. We find that JONSWAP data correspond to "semi-stable" 
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Figure 2: (a) Nonlinear spectrum and (b) evolution of Umax for JONSWAP data (7 = 4 and 
a = .016) that is near homoclinic data. Dashed curve corresponds to 2.2Hs. 



N-phase solutions, i.e. we interpret the data as perturbations of A^-phase solutions with one 
or more unstable modes (compare Fig. 2(a) with the spectrum of an unstable A^-phase 
solution in Fig. 1). In Fig. 2(a) the splitting distance 5(A+,A_) ~ .07, while in Fig. 3(a) 
5(A+,A_) ~ .2. Thus the JONSWAP data can be quite "near" homoclinic data as in Fig. 
2(a) or "far" from homoclinic data as in Fig. 3(a), depending on the values of the phases 
(pn in the initial data. For all the examined values of a and 7 we find that, when a and 
7 are fixed, as the phases in the JONSWAP data vary, the spectral distance 6 of typical 
JONSWAP data from homoclinic data varies. 

Most importantly, irrespective of the values of the JONSWAP parameters a and 7, in 
simulations of the NLS equation (0) we find that extreme waves develop for JONSWAP 
initial data that is "near" NLS homoclinic data, whereas the JONSWAP data that is "far" 
from NLS homoclinic data typically does not generate extreme waves. Figs. 2(b) and 3(b) 
show the corresponding evolution of the maximum surface elevation, Umax, obtained with 
the NLS equation. Umax is given by the solid curve and as a reference, 2.2Hs (the threshold 
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Figure 3: (a) Nonlinear spectrum and (b) evolution of Umax for JONSWAP data (7 = 4 and 
a — .016) that is far from homoclinic data. Dashed curve corresponds to 2.2Hs. 



for a rogue wave) is given by the dashed curve. Hs is the significant wave height and is 
calculated as four times the standard deviation of the wave amplitude. Figure 2(b) shows 
that when the nonlinear spectrum is near homoclinic data, Umax exceeds 2.2Hs (a rogue wave 
develops at about t — 40). Figure 3(b) shows that when the nonlinear spectrum is far from 
homoclinic data, Umax is significantly below 2.2Hs and a rogue wave does not develop. In 
this way, we correlate the occurrence of rogue waves characterized by JONSWAP spectrum 
with the proximity to homoclinic solutions of the NLS equation. 

The results of hundreds of simulations of the NLS equation consistently show that prox- 
imity to homoclinic data is a crucial indicator of rogue wave events. For example. Fig. 4 
shows the synthesis of 200 random simulations of the NLS equation for JONSWAP initial 
data for different (7,0;) pairs (with 7 = 2,4,6,8, and a — .012, .016). For each such pair 
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Figure 4: Strength of Umax/Hg vs. the splitting distance 5(A+, A_). 

(7, a), we performed 25 simulations, each with a different set of randomly generated phases. 
Each circle represents the strength of the maximum wave {Umax/Hs) attained during one 
simulation as a function of the splitting distance 5(A+, A_). The results for the particular 
pair (7 = 4, q; = 0.012) is represented with an asterisk. A horizontal line at Umax/Hs = 2.2 
indicates the reference strength for rogue wave formation. We identify two critical values 
61 = 0.08 and 62 = 0.22 that clearly show that (a) if 5 < Si (near homoclinic data) rogue 
waves will occur; (b) if 5i < 5 < 62, the likelihood of obtaining rogue waves decreases as S 
increases and, (c) if 5 > ^2 the likelihood of a rogue wave occurring is extremely small. 
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This behavior is robust. As a and 7 are varied, the strength of the maximum wave and the 
occurrence of rogue waves are well predicted by the proximity to homoclinic solutions. The 
individual plots of the strength vs. 8 for particular pairs (7, a) are qualitatively the same as 
in Fig. 4 as can be seen by the highlighted case (7 = 4, a = 0.012). These results give strong 
evidence of the relevance of homoclinic solutions of the NLS equation in investigating rogue 
wave phenomena for more realistic oceanic conditions and identifies the nonlinear spectral 
decomposition as a simple diagnostic tool for predicting the occurrence and strength of rogue 
waves. Finally we remark that the nonlinear spectral analysis can be implemented for other 
general data (non-JONSWAP) in order to predict the occurrence of rogue waves. 

This work was partially supported by NSF Grant No. ■nsf-dms0204714. 
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